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THEORETICAL  ANALYSIS  OF  DNA  LOW-FREQUENCY  PHONON  MODES 

M.  Bykhovskaia 

University  of  Virginia,  Dept.  Mol.  Physiology  &  Biol.  Physiscs 


Objectives: 

To  calculate  far  IR  absorption  spectra  of  DNA  double  helixes  with  known  nucleotide  sequences 
and  to  analize  the  affect  of  nucleotide  sequences  on  IR  absorption. 


Approach: 


DNA  potential  energy  is  calculated  and  minimized  using  a  semi-empiric  potential  for  interatomic 
interactions  (molecular  mechanics  method).  Normal  modes  are  calculated  as  eigenvectors  and 
eigenvalues  of  the  matrix  of  second  derivatives  of  molecular  energy.  The  ability  of  each  normal 
mode  to  add  to  the  DNA  absorption  spectra  is  evaluated  by  analyzing  of  charge  motions 
produced  by  the  corresponding  vibration.  The  dependence  of  the  absorption  coefficient  on  the 
wavelength  is  calculated  from  normal  vibrations  of  dipole  moments  of  the  molecules. 


Research  focus  and  impact 

Far-IR  spectroscopy  has  a  potential  to  become  a  detection  tool  for  biological  war-fare  agents. 
Nucleic  acids  are  bearers  of  genetic  information,  therefore  their  detection  is  important.  The  long¬ 
term  of  the  project  is  to  create  a  theoretical  basis  for  detection  of  certain  patterns  of  DNA 
structure  and  topology,  i.e.  to  detect  a  signature  of  particular  DNA  structures  in  its  far-IR 
absorption  spectrum.  The  present  project  is  intended  to  create  a  classification  of  DNA  low- 
frequency  vibrations  in  terms  of  their  correlation  with  nucleotide  composition,  three-dimensional 
structure,  and  ability  for  radiation  absorption.  Such  a  classification  is  fundamental  for  the 
understanding  of  DNA  biological  function,  because  the  function  of  DNA  crucially  depends  on  its 
dynamics.  In  addition,  detection  of  DNA  anomalies  is  important  for  the  understanding  of  genetic 
diseases.  As  the  first  step  for  this  long-term  goal,  we  calculated  and  analyzed  low-frequency 
vibrations  of  a  DNA  with  know  nucleotide  composition  and  far-IR  absorption  spectra. 


Technology  transfer  and  research  collaborations 

The  software  for  molecular  modeling  and  normal  mode  analysis  (program  packages  JUMNA  and 
LIGAND)  was  obtained  from  the  Lab  of  Dr.  Lavery,  Institute  of  Biochemistry  and  Physical 
Chemistry,  Paris,  France.  The  procedure  for  DNA  normal  mode  analysis  was  originally 
developed  by  Drs.  K.  Zacrzewska  and  Tap  Ha  Duong  (Dr.  Lavery’ s  Lab),  who  consulted  us  about 
the  usage  of  this  software. 
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The  method  for  evaluation  of  radiation  absorption  by  different  resonance  frequencies  was 
proposed  by  Dr.  B.Gelmont,  Dept,  of  Electrical  Engineering,  University  of  Virginia. 

The  farther  investigation  of  DNA  spectroscopy  is  planned  in  collaboration  with  Dr.  T.Globus, 
Dept,  of  Electrical  Engineering,  University  of  Virginia. 


Concepts  and  main  results 

We  calculated  vibration  frequencies  of  double  helical  DNA  fragment  (TA)6  and  and  (A)i2 
and  evaluated  optical  properties  of  all  thevibrations.  The  above  DNA  sequences  were  chosen 
because  its  far  IR  absorption  spectrum  have  already  been  measured  (range  cm'1  Powell  et  al 
1987). 

The  dependence  of  permittivity  8  on  radiation  frequency  to  can  be  described  by  the 
Kramers  -  Heisenberg  dielectric  function: 

e(<o)=6«  +  4rt  N  £  S  k/3(©k2-<a2  -  i  %  to/2  it  )  ,  (1) 

k 

where  index  k  denotes  normal  modes,  cok  are  vibrational  frequencies,  S  k  are  oscillator  strengths, 

Yk  are  oscillators  dissipations,  and  N  is  the  number  of  molecules. 

Absorption  coefficient  (a)  is  proportional  to  the  imaginary  part  of  permittivity: 

a=  ( to/  nc)  Im  (s ),  (2) 

hence  its  dependence  on  frequency  v  (v=<o/2ji)  is  determined  by: 

a(  v )~  v2  I  SkYk / ( (vkV) 2  +  (yk  v) 2)  (3) 

k 

Normal  mode  analysis  enables  calculations  of  normal  frequencies  v*  and  oscillator  strengths  S  k ,  while 
oscillator  dissipations  Yk  remain  to  be  determined  from  the  comparison  of  theory  and  experiment.  In  the 
first  approximation,  however,  the  decay  can  be  considered  frequency  independent.  In  this  simplified 
form  the  absorption  will  depend  on  the  frequency  as: 

a(  v  )  ~  v2  y  E  Sk  /  ( (vk2-v2) 2  +  (yv)  2)  (4) 

k 

The  realistic  estimation  of  y  obtained  from  the  measurements  of  (poly  dAXpoly  dT)  absorption  spectra  is  from  2  to 
7  cm  [10]. 


Standard  B-helical  conformations  were  created  by  the  program  package  JUMNA  [15-17]. 
Conformational  energy  was  calculated  as  a  sum  of  an  energy  of  van  der  Waals  and  electrostatic 
interactions,  torsion  rotation  potentials,  harmonic  potentials  of  bond  angle  deformations,  and  the 
energy  of  hydrogen  bonds  deformations  using  FLEX  force  filed  [18],  The  energy  was  minimized 
in  the  space  of  internal  coordinates  of  a  molecule  (torsion  and  bond  angles)  using  software  the 
package  LIGAND  [19]. 

Normal  modes  (eigenfrequencies  W  and  eigenvectors  A)  were  calculated  in  the  space  of  internal 
variables  as  solutions  of  the  equation: 
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HAW  =  FA 


(5) 


where  W  is  a  diagonal  matrix  with  elements  <Dk2  (vibrational  frequencies);  F  is  a  matrix  of  second 
derivatives  of  potential  energy  and  H  is  a  matrix  of  second  derivatives  of  kinetic  energy.  The 
normal  modes  (matrixes  A  and  W)  were  calculated  as  described  in  [71  using  the  program 
LIGAND. 

Each  normal  vibration  produces  fluctuations  of  the  dipole  moment  of  a  molecule.  Oscillator 
strength  Sk  corresponding  to  k  normal  mode  can  be  expressed  as  a  square  of  normalized  amplitude  (p) 
of  the  dipole  moment  deviation  [20]: 

Sk  =  (pk)2  (6) 

The  normalized  dipole  moment  p  can  be  expressed  as: 

p  =  I  ei  a*  /  Vmj  (7) 


where  index  i  denotes  atoms,  ej  are  partial  atomic  charges,  m;  are  atomic  masses,  and  aj  are  eigenvectors 
(atomic  displacements)  normalized  in  such  a  way  that: 

£(ai)2=l  (8) 

i 

The  equation  (7)  was  obtained  [20]  under  the  classical  mechanical  approximation  ( hv«kgT). 

However,  we  show  in  the  Appendix  that  (6)  is  also  applicable  in  ultraquantum  approximation 
(hv»kBT). 

The  procedure  consists  of : 

1)  Initial  optimization  of  the  B-helix  by  the  program  JUMNA; 

2)  Energy  minimization  and  calculation  of  normal  modes  (eigenfrequencies  cok  and  eigenvectors 
Ak )  by  the  program  LIGAND. 

3)  Calculation  of  three  dimensional  structures  of  a  molecule  along  normal  vibrations  [7]. 

4)  Calculation  of  oscillator  strengths.  The  variable  dipole  moment  of  a  molecule  (P)  is 
calculated  as: 

P=Sej  i-j  (9) 

Where  the  index  i  denotes  atoms  and  the  vector  r,  is  a  displacement  of  i*  atom  from  equilibrium. 
Oscillator  strengths  Sk  are  calculated  as: 

Sk  =  (Pk  )2max  /  2  mj  (n)2  ( i  o) 

where  Pk  is  produced  by  k01  vibration. 

5)  Calculation  of  absorption  spectra  a(  v )  (Eq.  4). 

Oligonucleotides  (Poly  dA)12(Poly  dT)12  and  (Poly  dAdT)6  (Poly  dTdA)6have  360 
degrees  of  freedom  in  the  space  of  internal  coordinates  of  a  molecule  (torsion  and  bond  angles). 
Respectively,  360  normal  modes  were  found  for  each  sequence.  Their  frequencies  lie  below  900 
cm-1  (Fig.  1),  so  there  is  almost  no  overlap  with  vibrations  of  covalent  bonds  which  have 
frequencies  above  750  cm  1  [21].  The  two  oligonucleotides  have  similar  mode  densities  in  the 
region  above  250  cm'1 ,  which  essentially  involve  vibrations  of  valence  angles  [7].  It  could  be 
expected  that  the  modes  that  reflect  vibrations  of  torsion  angles  would  be  more  sensitive  to  the 
conformation  and  flexibility  of  the  helix.  Indeed,  the  calculated  density  spectra  of  the  two  oligos 
differ  markedly  in  the  region  below  250  cm'1.  Specifically,  (Poly  dAdT)6  (Poly  dTdA)6  has  very 
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distinct  peaks  around  15  and  175  cm'1 ,  while  homoplymer  (Poly  dAXPoly  dT)  with  unique 
conformational  properties  [14]  has  a  broad  spectrum  with  no  pronounced  peaks  (Fig.  1).  The 
earlier  normal  mode  analysis  of  (Poly  dAdT)$  (Poly  dTdA)$  performed  with  rigid  valence  angles 
[7]  revealed  only  one  peak  at  15  cm’1 ,  which  indicates  that  the  second  peak  at  175  cm*1 
essentially  involves  vibrations  of  valence  angles. 
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IQ¬ 
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Poly  (dAdT)  Poly  (dtdA) 


Vibrational  frequency  (cm'1) 


Fig.  1.  Density  of 
normal  modes. 
Histograms  on  the  left 
represent  all  the 
calculated  modes.  The 
scale  is  expanded  (right 
panels)  to  show  the 
lowest  frequency 
modes. 


Both  Poly( 
dA)Poly(dT)  and 
Poly(dAdT)Poly(dTd 
A)  have  a  resonance 
absorption  at  the 
frequency  60  cm'1 
[10],  therefore  we 
examined  deviations 
of  a  dipole  moment 
Pk  (Eq.  7)  produced 
by  the  strongest  mode 
with  the  frequency 
near  60  cm*1  (Fig.  2). 
This  mode  was 


previously  assigned 

to  vibrations  along  the  helix  axis  (Z),  however,  Fig.  2  demonstrates  that  the  traverse  component 
Px  has  comparable  amplitude. 


Fig.  2.  Fluctuations  of  the  dipole  moment  Pk  of  the  molecule 
(Poly  dAdT)6  (Poly  dTdA)6  corresponding  to  a  normal 
vibration  with  a  frequency  60.3  cm'1  at  the  room 
temperature.  The  double  helix  is  oriented  along  the  axe  Z. 
The  amplitude  of  (P*)2  fluctuation  is  used  in  Eq.  10. 

The  amplitude  of  P2  fluctuations  determines  the 
strength  of  each  vibration  (Eq.  9).  Normal  modes  of 
the  two  examined  oligos  demonstrated  very  different 
distributions  of  oscillator  strengths.  As  would  be 
expected,  for  the  molecule  with  homogeneous  strands 
Poly(dA)Poly(dT)  the  oscillator  strengths  appeared  to 
be  very  uniform,  while  heterogenity  in  the  strands 
composition  introduces  many-fold  variance  (Fig.  3). 
The  two  lowest-frequency  modes  (2.1  and  2.6  cm'1) 
have  strengths  3-4  times  greater  than  the  average 
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Absorption 


strength  for  the  other  modes.  These  two  lowest-frequency  vibrations  are  have  strong  P 
components  perpendicular  to  the  helical  axe,  and  they  apparently  reflect  helical  bending  [7, 22]. 
The  strongest  modes  of  the  heteropolymer  (Fig.  3  A)  also  lie  in  the  lowest  frequency  range 
(mostly  below  50  cm'1).  This  result  is  not  surprising,  for  the  strongest  modes  reflect  the  relative 
motions  of  the  strands,  which  in  turn  involve  low-frequency  torsion  rotations.  Correspondingly, 
the  normal  modes  spectrum  weighted  by  oscillator  strengths  is  essentially  shifted  towards  low- 
frequencies  (Fig.  4A  compared  to  Fig  2A)  for  the  heteropolymer,  while  for  the  homopolymer 
heterogenity  in  oscillator  strengths  has  a  minor  effect  compared  to  the  normal  mode  density 
itself. 


Frequency  (cm*1) 

Fig.3.  Distributions  of  oscillator  strengths  for  normal  modes.  Fig.  4.  Density  of  normal  modes  weighted 

with  their  oscillator  strengths.  Data  are 
normalized  to  represent  the  oscillator 
strength  per  normal  mode 


We  calculated  absorption  spectra  of  the  two  molecules  (Fig.5)  according  to  Eq.  4  using 
two  reasonable  approximations  for  the  band  width  y  obtained  experimentally  (2  and  7  cm-1, 
[10]).  The  positions  of  resonance  peaks  do  not  depend  on  the  band  width  (Fig.  5),  and  in  the 
longer  wave  range  (below  220  cm  *)  they  are  remarkably  similar  to  those  observed 
experimentally  (Table  1). 


Poly(dA)Poly(dT)  Poly(dAdT)Poly(dTdA) 

y=2  cm  . 


Fig.  5.  Absorption  spectra  (arbitrary  units) 
of  Poly(dA)12Poly(dT)i2  (left)  and  (Poly 
dAdT)e  (Poly  dTdA)6  (right)  calculated  at 
two  values  (2  and  7  cnr)  of  the  band 
width  y. 
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Table  1.  Resonance  frequencies  for  Poly(dA)Poly(dT).  All  numbers  are  in  units  cm'1 


Calculations _ Experiment  [10] 


20 

unknown 

43 

unknown 

57 

62 

71 

80 

94 

95 

106 

106 

131 

136 

160 

- 

172 

170 

199 

- 

211 

214 

- 

238 

270 

~ 

348 

. 

400 

- 

459 

- 

Another  similarity  of  the  calculated  and  experimental  spectra  is  the  steep  fall  in 
absorption  intensities  as  the  frequency  exceeds  250  cm'1  for  both  sequences.  This  feature  is 
caused  by  the  decrease  in  normal  mode  densities  and  is  not  related  to  their  optical  properties. 
Finally,  Poly(dAdT)Poly(dTdA)  demonstrates  higher  absorption  at  lower  frequencies  (below  100 
cm' )  due  to  the  higher  oscillator  strengths  of  its  lower  frequency  modes,  while 
Poly(dA)Poly(dT)  has  a  plateau  between  20  and  200  cm'1 .  This  difference  in  the  spectra  has 
also  been  observed  experimentally  [10]. 

Our  results  indicate  a  strong  resonance  peak  around  20  cm  1  for  the  both  sequences. 
While,  no  experimental  data  on  IR  absorption  is  available  in  this  frequence  range,  Raman 
scattering  [24]  clearly  indicates  resonance  peaks  around  20  cm'1  (16.2,  193.  3  and  23.3  cm'1) 
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